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We obtain exact solutions of the (2+1) dimensional Dirac oscillator in a homogeneous magnetic 
field within the Anti-Snyder modified uncertainty relation characterized by a momentum cut-off 
(p < p m ax = 1/y/P). In ordinary quantum mechanics (/3 —> 0) this system is known to have a single 
left-right chiral quantum phase transition (QPT). We show that a finite momentum cut-off modifies 
the spectrum introducing additional quantum phase transitions. It is also shown that the presence 
of momentum cut-off modifies the degeneracy of the states. 
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I. INTRODUCTION 

In recent years quantum mechanical models based on 
modified commutation relations (MCR) have been stud¬ 
ied by many authors. Usually MCR’s involve a minimum 
position uncertainty leading to a minimal length mm 
( Hy/j3 ). But there are also MCR’s which lead to a maxi¬ 
mum momentum [SHIP] (l/i//3). A number of quantum 
mechanical models in the presence of a maximum mo¬ 
mentum parameter have also been investigated mm- 
The main reason for studying quantum mechanical mod¬ 
els with MCR is that the spectrum gets modified because 
of the presence of a minimum length/maximum momen¬ 
tum parameter which may eventually be detected in low 
energy experiments. 

The Dirac oscillator DEH23 in the presence of a homo¬ 
geneous magnetic field is one of the few exactly solvable 
problems in relativistic quantum mechanics |22j . This 
combined system is known to exhibit a quantum chiral¬ 
ity phase transition whenever the magnetic field strength 
exceeds a critical value [33H35] . Such system has received 
quite a lot of attention also in view of the possible appli¬ 
cations to the physics of recently discovered materials like 
graphene [25H25] . germanene (991 [3(1] and silicene [5TH51] . 
Indeed it has been known that, in these materials, at the 
K, K' points of the Brillouin zone the charge carriers are 
described by an effective 2 dimensional Dirac equation. 
The Dirac oscillator coupling has been proposed [35. to 
arise from the interaction of the charge carriers with an 
effective internal magnetic field due to collective motion 
of the electrons in the planar hexagonal lattice of the 
carbon atoms. 

This system has also been studied in the presence of a 
minimal length |36l as well as in non-commutative space 
[371 [Ml- In all these cases it was shown that there are a 
multitude of quantum phase transitions depending on the 
minimal length/non-commutative parameters. In this 
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work our objective is to study again the same system 
but in the Anti-Snyder space i.e, under a MCR which 
incorporates a momentum cut-off. One reason why we 
have taken up this model is that in contrast to the mini¬ 
mum length scenario in the present case the presence of 
a maximum momentum parameter changes the nature of 
the spectrum from an infinite one to a finite one. Con¬ 
sequently it is of interest to examine how the quantum 
phase transitions are influenced by the maximum mo¬ 
mentum parameter. 

We will show that in the Anti-Snyder model the (2+1)- 
dimensional Dirac oscillator of frequency to in an external 
magnetic field (_B 0 ) admits a finite and bounded spectrum 
(of bound states). Similarly to what has been found in 
ref. [35] we find that the excited energy levels can be 
classified in two classes one of which disappears in the 
ordinary quantum mechanical limit, /3 -+ 0, and whose 
corresponding states are characterised by the presence of 
quantum phase transitions which accumulate towards the 
critical field Bq = 2ul ^ Ic . The other class of levels does 
not disappear in the ordinary quantum mechanical limit, 
/3 -A 0, and indeed reduces to the known spectrum of the 
system in this limit. We postpone the discussion of the 
continuum states to a future work. 


We would also like to recall that recently the one di¬ 
mensional version of the Dirac oscillator has been realised 
in the laboratory [33]. In addition practicable prospects 
of realising soon the two dimensional version of the Dirac 
oscillator have also been reported in the literature [39b 
l4l] . The exact solutions presented here may have then 
a direct relevance to the physics of experimentally acces¬ 
sible 2-dimensional systems which have received a lot of 
attention recently in the literature like [261134] . 


The rest of the paper is organised as follows: Sec. [TT] 
introduces the Anti-Snyder model and defines the prob¬ 
lem; Sec. |III| describes the transformation of the problem 
to a Poschl-Teller potential; Sec. m and Sec. [V] present 
the explicit solutions respectively in the case of weak and 
strong magnetic field. Sec. VI gives a discussion of the 
main properties of the spectrum and the quantum phase 
transitions. Finally Sec. m presents the conclusions. 
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II. FORMULATION OF THE PROBLEM 

Let us first present the basic facts about the Anti- 
Snyder space. The commutation relations between the 
coordinates (a;*) and momenta {p 3 ) are given by [IT]: 

[xi,Xj\ = 2ih/3(l - /3p 2 )e ijk L k , \pi,Pj] = 0, 

[Xi,pj] = ihSiji 1 - dp 2 ), ]pu L 3 \ = iheij k p k , 

, Lj\ — [L ? ;. Lp — ihCijfcLfc, 

(1) 

where Li = 1 _ 1 /3 p 2 CijkXjPk • A representation of the oper¬ 
ators Xi and pi can be taken as 


d 

Xi = ih(l - /3p 2 ) — , pi=Pi 

dPi 


( 2 ) 


It may be pointed out that in the present case the scalar 
product is defined as 


III. TRANSFORMATION TO A 
POSCHL-TELLER II POTENTIAL 

In order to solve the eigenvalue problem in Eq. © we 
proceed as follows: 

V’l 1 ) = ( p ), <A (2) = e« m+1 *V 2) (p) ( 10 ) 

where m = 0, ±1, ±2,... denotes the angular momentum 
quantum number in the p space. 

Now applying the operator H from the left on Eq. ([5]), 
the eigenvalue equations for the components can 

be found as: 


P_P+V> (1) = e 2 ^ (1) , 
P+P_v> (2) = e 2 ?/> (2) , 


e 2 = 


E 2 - M 2 c 4 e + e_ 


(11a) 

(lib) 

(11c) 


( f\g) = 




d 2 p 


(1 - (3p 2 ) 


f*(p)g(p ) 


( 3 ) 


Let us now consider the Hamiltonian corresponding to 
(2 + 1) dimensional Dirac oscillator in the presence of a 
homogeneous magnetic field B = (0,0, Po)- In the Anti- 
Snyder space the corresponding Hamiltonian is given by 

H = ccr ■ (p — iMua z x + - A) + a,Me 2 , (4) 

c 

where er = (a x ,a y ) and A = (—B 0 y/2,B o x/2,0) denote 
the Pauli matrices and the vector potential respectively. 
We shall now solve the eigenvalue equation corresponding 
to the Hamiltonian Q. The eigenvalue problem reads 


Hp = 


M c 2 cP 
cP+ — Me 2 


^ (1) \ - E 
pW )~ E 


( V’ (1) \ 

V v> (2) ) 


or 


cP-ip^ = e_ 
cP+'ip^ = e + 


( 5 ) 


(6a) 

(6b) 


where e± = E ± Me 2 and the operators P± are given by 

P+ = e+" 

P_ = e~ m 


p-\(l-(3p 2 ) (dp + ^do, 


P + A (1 - dp 2 ) [d p - 


( 7 ) 


with 


A = — ° — Mhui = Mh (oj c — u>) , 
2c 


eB 0 
2 Me 


( 8 ) 


p x =p costf, Py=p sind, p 2 x +p 2 y =p 2 (9) 

From Eq. ^ it is seen that depending on the strength 
of the magnetic field (in comparison to the oscillator 
strength), A can either be positive or negative. 


With a little algebra then the above Eqs. (11 a|TIbt can 
be put in the following form: 

- f(p) d f 2 +9(p) d ^j— + /ii(p)v? (1) = (e 2 //?A 2 V (1) 

dp z dp 

( 12 ) 

- fip) d ? 2 + 9 , (p)-^— + h 2 {p)p (2) = (e 2 //? \ 2 )p {2) 

dp z dp 

(13) 

where 


fip) = 


(1 - dp 2 f 


g(p) = (! - dp 2 ) 


2 P- 


(1 - dp 2 ) 
dp 


, p 2 (1 — (3p 2 )(2m + 2 — 2m/3X) 

hM = W + - W~ 


(1 — f3p 2 ) 2 m 2 


j3 p 2 

v 2 (1 — 0v 2 ) 

h 2 {p) = + (2m + 2m/3A + 2/3A) -——- 

(m + 1) 2 (1 — dp 2 ) 2 
/3p 2 


(14) 


(15) 


In order to solve the equations (12) and (13) we per¬ 
form the transformations 

<P (1,2) (P) = P(?#1, 2 (P), 9 = J^=/P ( 16 ) 


Pip) = e /x(p) dp , X(P) = ^^ (17) 

namely 

¥> (1,2) (p) =P _1/2 0i,2(p), 9 = tanlW^y^p) (18) 
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and obtain from Eqs.(12) and (13) 


r^ +m 


= k 2 (j>i(q), * = 1,2 


k 2 = 


l -VP 

/3X 2 


(19a) 

(19b) 


Eq. (19a) is effectively an one dimensional Schrodinger 


equation with a potential given by 

T _ , , (—2 + (3 A + 2ro/3A)(— 2 + 3/3A + 2m(3X) 2 

V i(®) =-^2- sech 9 

+ (m — ^)(ra + cosech 2 g, 0 < q < oo (20) 

r , v ( — 2 — /3A + 2mf3\)(—2 + /3A + 2mf3X) 2 

W =- ^2 - sech q 

1 3 

+ (m + -)(m + cosech 2 g, 0 < g < oo (21) 


It may be noted that the potentials (20) and (21) can be 


identified with the Poschl-Teller potential II of the form 


U(x) = —A(A + 1) sech 2 cc + B(B — 1) cosech 2 x, 
0 < x < oo, A > B > 0 


( 22 ) 


The above potential is exactly solvable. The eigenvalues 
and the corresponding (unnormalised) eigenfunctions of 
the Schrodinger eigenvalue problem as in Eq. ((19a)) are 
given by HS; 23]: 


Aj > Bj > 0 for j = 1,2 whose resulting constraints 
are shown in the third and fourth columns of Table [Hand 
Table HO 

Before discussing the explicit general solutions for the 
eigenvalues and wavefunctions based on the above results 
we conclude this section by a separate discussion of the 
zero modes, e.g. the states with energy E = ±Mc 2 (or 
E = 0 in the massless case). Such states can be identi¬ 
fied solving explicitly the first order equations in Eq. © 
which read: 


cP_^ (2) = (E- Me 2 ) ip {1) (24a) 

cP+ip {l) = (E + Me 2 ) ip {2) (24b) 


From Eqs. (24) it is seen that the eigenfunctions cor¬ 
responding to the eigenvalues E = ±Mc 2 are correctly 
indentified by the following set of conditions: 


E = +Mc 2 , -a ip {2) = 0, P + V> (1) = 0 (25a) 

E = -Me 2 , -a ^ (1) = 0, P-ip {2) = 0 (25b) 


The corresponding first order differential equations are 
readily solved explicitly. Using Eq. ([3| it can be s hown 
that the acceptable (normalizable) solutions of Eq. (25a) 
are given by 


ip = Ne 


imO ( ty? 


o 




l 1 ) = p m (1 - pjp 2 ) 


E = +Mc 2 , 
m > 0, A < 0, (26) 


k 


2 

n 


(pn 


-{A - B -2n) 2 , 


n = 0,1, 2,... < 


A-B 

2 


/V(cosh2x — l) B / 2 (cosh2:r + 1) A ^ 2 x 
2 Pi(-n, B - A + n; B + 1/2; 1 ~ co 2 sh2a: ) 


(23a) 

(23b) 


where 2 Pl(— n, B — A+n; P + 1/2; 


1—cosh 2x 


) is the hyper¬ 


geometric function, and N a normalisation constant. 
Note that the conditions A > B > 0 are required for 
the wave functions to be normalizable. Now identifying 
the potential (21) with the one in (22) we have to solve 
the following two systems of quadratic equations in the 
unknown quantities (Ai,Bi) and (A 2 ,B 2 ): 


(—2 + /JA + 2m/3A)(—2 + 3/3A + 2mj3\) 

+1} = - Wx 2 -’ 

Pi(Pi - 1) = (m - i)(m + i); 

^ (—2 — /3A + 2m/3A)(—2 + /3A + 2mf3\) 

A 2 (A 2 + 1) --^2y2-, 

B 2 (B 2 - 1) = (m + *)(m + ?). 

Each system admits four solutions (a, 6, c, d) which are 
shown in Table [I] and Table |TT] (first and second columns). 
In order to ensure that the resulting wave functions are 
normalizable it is necessary to impose the conditions 


where N is a normalization constant. 
Similarly Eq. (25b) is solved by: 


ip = /Ve*( m+1 ) 1? ^ ^(2) ^ , E = —Me 2 , 

=p-( m+1 \i - /3 p 2 )^h m < —1, A > 0, 


(27) 


where N is again a normalization constant. We then 
realize that when A < 0 there is a zero mode at E = 
+Mc 2 while for A > 0 the zero mode is in the negative 
branch of the spectrum. In both cases the zero mode is 
infinitely degenerate with respect to the angular quantum 
number in. 


IV. WEAK MAGNETIC FIELD (w c < w, A < 0) 

Having discussed the zero modes explicitly in the 
previous section we now turn our attention to the general 
excited states, starting from the region of negative values 
of A. 


(*) States of finite degeneracy. 

We start considering solutions (a) in Table I] and Ta¬ 
ble |nj The reduced eigenvalues k 2 (last column), for 
both components </>i j2 are obtained using the solutions 
of Ai,Bi, and A 2 ,B 2 (third and fourth columns) in the 
eigenvalue relation in Eq. (|23a). In this case we have 
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Ai 


Bi 

Ai > Bi 

Bi > 0 

A < 0 

0 < A < yg 

rv n,m 

(a) 

1 

/3A 

+ 1 

+ TO 

-TO+ I 

m >Wx 

TO < 0 

m <™<° 

- 

- A (~Wx +m-n) 

(b) 

1 

13 A 

+ 1 

+ TO 

m+\ 

-jx >0 

TO > 0 

TO > 0 

- 

- 4 (-2k~ n ) 

(c) 

1 

P A ' 

3 

2 

- m 

-TO+ | 

jx >2 

TO < 0 

- 

TO < 0 

~ 4 ( 2/3A — 4 ~ n ) 

(d) 

1 

/3A ' 

3 

2 

- m 

TO + | 

m <Wx~ 1 

TO > 0 

- 

°< m < m- 1 

- 4 {wx~ 1 - m - n ) 


Table I. We give here the possible solutions (a, b, c, d) of the coefficients A± and B\ in terms of the parameters m and j3 A when 
identifying the potential V\ (q) in Eq. [51] with the Poschl-Teller potential of type II, U{x ) in Eq. |22| We obtain the m ranges 
respectively for the cases A < 0 (fifth column) and A > 0 (sixth column) combining the results form the conditions Ai > B± 
(third column) and B i > 0 (fourth column). Cases (c) and (d) when A > 0 both imply the condition A < ^g. 




A 2 


b 2 

-A 2 > B 2 

B 2 > 0 

-h < A <o 

A > 0 

k 2 

n 'n,m 

(a) 

1 

/3A 

1 

2 

+ TO 

-TO - | 

m> Wx 

TO < 

-1 

Wx <m <~ 1 

- 

- 4 (-^a +m-n) 

(b) 

1 

/3A 

1 

2 

+ TO. 

TO+ I 

JX<- 2 

TO > 

-1 

TO > — 1 

- 

~ 4 (~ 2^A — 4 — n ) 

(c) 

1 

0 A " 

1 

2 

- m 

-TO - 1 

M >0 

TO < 

-1 

- 

TO < —1 

~ A (wx ~ n ) 

(d) 

1 

13 A " 

1 

2 

- m 

TO+f 

m <m~ 1 

TO > 

-1 

- 

—! < TO < 2 ^ — 1 

~ 4 {wx -i-TO-n) 


Table II. The same as in Table [I] but for the potential V 2 ( 9 ) in Eq. 20 Cases (a) and ( 6 ) when A < 0 both imply the 
condition — ^ < A. Note that when combining the upper and lower solutions the possible values of A are limited to the interval 
— < A < 2 ^ 3 . This can be qualitatively understood since A defined in Eq. ^1 has the dimension of a squared momentum. 

Hence, being the Anti-Snyder model characterised by a maximum momentum p m ax = it follows a bounded allowed range 

in A. 


A 1 — Bi = A 2 — B 2 = —1/(2/5A) + to so that the eigen¬ 
values of the two components match using for each one 
of them the same radial quantum number n. We must 
however take the intersection of the ranges of the angular 
quantum number to from Table [I] and Table [n| In this 
context we note that 

{ 2^A < m — 0 } n { < to < — i} 


{2 


< TO < —1 


Note also that this implies in particular < — 1 which, 
since here A is negative, translates into a global condi¬ 
tion — jp < A (i.e. valid for both components). All neg¬ 
ative values A < — W which would have been allowed for 
the solution of the upper component (Table [I]) must be 
excluded when matching with the solution of the lower 
component which requires — js < A. 

We then finally get: 


k l,m = +m-nj 

n = 1, • 


< - 


2/3A 


TO, 


i 2/3 < 0 

m <m - _1 


(28) 


with the eigenvalue E n m of the Dirac problem computed 
via Eqs. ( 23a|19b ) and the corresponding upper and 


lower components of the spinor solution with the same 
radial quantum number n. 

(ii) States of infinite degeneracy. 

We now turn our focus on solutions ( b ) of Table [T] and 
Table [TT] The reduced eigenvalues k ^ (last column), for 
both components (f>i t2 are obtained again using the solu¬ 
tions of and A 2l B 2 (third and fourth columns) 

in the eigenvalue relation in Eq. |23a[ Now we find 
Ai — Bi = —-gh and A 2 — B 2 = — 2. We see 

then that the eigenvalues (last column) do not match if 
we use for both components the same radial quantum 
number n. This implies that in writing a spinor solu¬ 
tion of the Dirac problem we must connect components 
with the radial quantum numbers shifted by one unit: 
e.g. with ^ 2 2i_ In addition we must take as accept¬ 
able only the intersection of the ranges of the angular 
quantum number to from Table [T] and Table [TT} Thus: 

{to > 0} n {to > — 1} => {to > 0} 

As in the case of the solutions (a) while the condition 
in the third column of Table [I] is always satisfied (being 
A < 0) the same condition on Table [TTj < — 1, implies 

— ^ < A translates into a global condition (i.e. valid 
for both components). Again, all negative values A < 

— Ap which would have been allowed for the solution of 
the upper component (Table [T| must be excluded when 
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matching with the solution of the lower component which 
requires — A < A. 

We then finally get: 

kn,m ~ —4 (~ 2/3A _ n ) ’ ~ 2/3 < ^ ^ (29) 

n = 1, • • • < — 2 ^ + m, m> 0 
with the eigenvalue E n , m of the Dirac problem computed 


via Eqs. ( 23a|19b ) and the corresponding upper and lower 
components of the spinor solution with the radial quan¬ 
tum number n shifted by one unit, e.g. with 

The full spinor solutions in the region — ^ < A < 0, 
corresponding to cases (a) (finite degeneracy, Eq. (28)) 


and ( b ) (infinite degeneracy, Eq. (29)) illustrated in the 
above detailed discussion, are explicitly written down in 
the following Eqs. ( |30l3l| ), in terms of unspecified nor¬ 
malisation constants C_, C__. 


E n ,m = CyJ M 2 c 2 + | - 4/3A 2 (-^y + m- n) 2 , tp n>m = C_ ^ 

<ph]m =p Bj ~ 2(1 - ftp 2 ) 2 2 Fi(— n, Bj - Aj +n-,Bj + 1/2; — Hl> ), sgr < m < -1, j = 1,2 


e + +imi? (1) 

2mA 

e +i(m+i)^(2) m 

pp 2 


, n 1,.... < 


2/3A 


1 — (3p 2 ’ 2/3A 


En.m = cy / M 2 C 2 + g - 4/3A 2 (— - «) 2 , V’n.m = C-_ 


+imt) (1) 

£ V^n,m 


0 +i(m+l)i9 


e_ 


. .. 1 A j B j Q ny2 

Vn,m = p Bj ~ 2 (1 - /V) 2 2 Fi(-n, - A, + n; + 1/2; ~ 1 _ 2 ), m >0, j = 1, 2 


2(m + 1)A 

Pp 


( 2 ) 

V’n-l.n. 


(30) 

, n = 1, .... < — 2 / 3 A 

(31) 


V. STRONG MAGNETIC FIELD (w c > w, A > 0) 

We now discuss the spectrum in the case when the 
magnetic field is such that the corresponding cyclotron 
frequency uj c is larger than the Dirac oscillator frequency 
ui which amounts to requiring A > 0. 

While the procedure is quite similar to the one exposed 
in the previous section the corresponding spectrum turns 
out to have important differences relative to the previous 
case. Such differences are responsible for the quantum 
phase transition(s). 

(*) States of infinite degeneracy. 

We now turn our focus on solutions (c) of Table [T] and 
Table |h| The reduced eigenvalues k 2 (last column), for 
both components are obtained again using the solu¬ 
tions of and A 2l B 2 (third and fourth columns) 

in the eigenvalue relation in Eq. |23a| Now we have 
Ai — Bi = — 2 and A 2 — B 2 = l/(/3A). We see then 

that the eigenvalues (last column) do not match if we 
use for both components the same radial quantum num¬ 
ber n. This implies that in writing a spinor solution of 
the Dirac problem we must connect components with the 
radial quantum numbers shifted by one unit: e.g. 

(2) 

with ipn '. Again we must take as acceptable only the in¬ 
tersection of the ranges of the angular quantum number 
to from Table HJ and Table HT1 Thus: 

{to < 0} n {to < —1} =$■ {to < —1} 


As in the case of the solutions (b) while the condition in 
the third column of Table[IT]is always satisfied (being A > 
0) the same condition on Table |T| > 2, implies A < ^ 

and translates into a global condition (i.e. which must 
be valid for both components). Again, all positive values 
A > pp which would have been allowed for the solution of 
the lower component (Table |n| must be excluded when 
matching with the solution of the upper component which 
requires A < 

We then finally get: 


— n 


n ~ 1, < 2/3A > 


— OO < TO < —1 


(32) 


with the eigenvalue E njm of the Dirac problem computed 
via Eqs. ( 23a|19b ) and the corresponding upper and 
lower components of the spinor solution with the radial 
quantum number n shifted by one unit, e.g. with 

( 2 ) 

I Pu¬ 


pil) States of finite degeneracy. 

We are now left with the solutions (d) in Table HI and 
Table |h| The reduced eigenvalues k 2 (last column), for 
both components (j)\ t2 are obtained using the solutions 
of Ai,Bi, and A 2 ,B 2 (third and fourth columns) in the 
eigenvalue relation in Eq. (|23a). In this case we have 
Ai — B ± = A 2 - B 2 = l/(2/3A) -2-2 m so that the 
eigenvalues of the two components match using for each 
one of them the same radial quantum number n. We 
must however take the intersection of the ranges of the 
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angular quantum number to from Table [T] and Table [TT| 
In this context we note that 


0:S " 1< 4“ 1 } n {“ 1Sra< 4~ 1 


We then finally get: 


k 2 

n 'n,m 


= 4 ( 2 ^a - 1 -m-n) 2 


n = i’— < m ~ l ~ m 


0<m< 2^-! 


(33) 




Note also that this implies in particular > 1 which, 
since here A is positive, translates into a global condition 
A < jp (i.e. valid for both components). All positive 
values A > ^ which would have been allowed for the 
solution of the lower component (Table [IT]) must be ex¬ 
cluded when matching with the solution of the upper 
component which requires A < 4. 


with the eigenvalue i? n ,m of the Dirac problem computed 
via Eqs. ( 23a|19b ) and the corresponding upper and lower 
components of the spinor solution with the same radial 
quantum number n. 

The full spinor solutions in the region 0 < A < 4 cor¬ 
responding to cases (c) (infinite degeneracy, Eq. (321) and 
(d) (finite degeneracy, Eq. (331) illustrated in the above 


detailed discussion, are explicitly written down in the fol¬ 
lowing Eqs. ( 34|35 ) in terms of unspecified normalisation 
constants C+, 


^4- 


En : m = C^M 2 C 2 + J - 4/?A 2 ( jgy - n) 2 , l />„, m = C+ 


Vn}m=p Bi 2(1 - Pp 2 ) 2 2 F 1 (-n, Bj - Aj + n; Bj + 1/2; - ^ 2 ), -oo < m < -1, j = 1,2 


Aj—Bj 


r +im'd £ + ,„(!) 

6 2mX^ n ~ 1,r ‘ 

+i(m+ 1)t? ..(2) 


Pp 2 


n=l 1 .... < 


/ 3p 2 


(34) 


En,m = CyjM 2 C 2 + | - 4/3A 2 (2^y 1 TO Tl) 2 , V’n.m = C+ 


e_ 


^m=p Bj 2(1 - Pp 2 ) 2 2 ^1 (“«, Bj - Aj + n\ Bj + 1/2; - ) ’ 0 < 




An —Bn 


2 (to + 1)A 

Pp 2 


p +imi9 (1) 

+z(m+l)tf (2) 

° rn,m 


n •■■■ < 2/3A 1 


ftp 2 


m < m ~ lj J = 1 > 2 


(35) 


VI. SPECTRUM AND QUANTUM PHASE 
TRANSITIONS: DISCUSSION 


It is well known that even in ordinary quantum me¬ 
chanics (/3 —> 0 - no maximum momentum -) when 
Cj c = oj there is a quantum phase transition. Here our 
main object is to study the impact of the Anti-Snyder 
model P ^ 0 with respect to the nature (and number) of 
such phase transition(s). 

Let us consider for instance the spectrum E nm as given 
in Eqs. ( 30|31 ) and Eqs. ( 34|35 ). Upon introducing the 
dimensionless variable: 


P = 


h(u> c — w) 
Me 2 


A 


M 2 c 2 ’ 


and the parameter: 


P = 


(36) 


(37) 


PM 2 c 2 ’ 

the energy eigenvalues from Eqs. ( |30pl| ) for A < 0, and 


in Eqs. ( 34|35 ) for A > 0 are easily obtained as: 


Eat = < 


±Mc 2 y 1 — 4 Up Af 

±Mc 2 J 1 + 4A fp 


p < 0 


p > 0 


(38) 


where for p < 0: A f = n + \m\ (negative to), Eq. (301, 
and A f = n, Eq. ([3l|); and for p > 0: A f = n, Eq. (341, 


and AT = n + m + 1 (positive to), Eq. (35). 

The energy eigenvalues E^j- are plotted as function of 
the p variable in Fig. [I] where it is quite evident that 
the discrete spectrum is indeed bounded: for each value 
of p there is a finite number of discrete energy levels and 
there is a maximal energy corresponding to the maximum 
momentum of the Anti-Snyder model. 

One first thing to notice about the spectrum is that it 
is always finite for any finite value of p and it is always 
bounded from above E < E max with E max easily com¬ 
puted from Eq. 38 as E max = Mc 2 y/1 + p* and straight 


forwardly related via Eq. m to the maximum momen¬ 
tum: E max = M 2 c 4 + j)^ OJ c 2 , as would be naturally 
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2 3 4 5 • • • • 0 ■■■•54 3 2 


Figure 1. (Color online) Positive branch of the discrete spectrum. The (positive) energy eigenvalues Em from Eq. (381, in 
units of Me 2 , are plotted as a function of the dimensionless variable p which is directly related to the intensity of the external 
magnetic field B. Note that for a given value of the dimensionless variable p we order the eigenvalues according to increasing 
values of the energy Em resulting in the (excited) energy levels E~£ which are assigned a color code according to the legend. 


expected for a relativistic particle with a maximum mo¬ 
mentum. 


For a given value of the dimensionless variable p we or¬ 
der the eigenvalues according to increasing values of the 
energy Em resulting in the (excited) energy levels E£ 
which are assigned a color code according to the legend, 
see Fig. [I] Since the functions Em{p) from Eq. (38), 
are non monotonic the level suffix k in general does 
not coincide with the suffix J\f which denotes instead 
a given function Em{p)- The energy levels are clearly 
separated in two distinct classes. One class of levels is 
built up (say for p > 0) with the decreasing branches of 
the Em(p) functions (k = 1,3,5,...). This levels, which 
have the peculiar property that they disappear in the 
[3 —> 0 limit ( p* -+ oo), suffer a discontinuity at each 


value 


p * 


-, j = 1,2 • • ■ which therefore denote a series 


2k+l+j ' 

of quantum phase transitions accumulating towards the 
p = 0 critical point. The class of levels (for p > 0) built 
out of the increasing branches of the functions Em{p) 
(ground state, k = 2,4,6 • • •) reduce in the /3 —>- 0 limit 
(p* —> oo) to the known standard result of ordinary quan¬ 
tum mechanics m extending to p —> oo. 


Therefore we can give the following definition of the 


two classes of excited levels: 


C 


< 


c 


> 


fc = 2, 4,6, • • • 

, (p < o) 

k = 1,3,5, • • • 

, (P > 0) 

fe = 1,3,5,--- 

, (P < 0) 

fc = 2, 4, 6, • • • 

i (p > o) 


(39) 


(40) 


We emphasise that in the vanishing /3 limit the class 
C> reduces to the ordinary spectrum of the (2+1)- 
dimensional Dirac oscillator while the class C< contains 
new levels relative to the ordinary spectrum which dis¬ 
appear in the limit ft -3- 0 (p* -+ oo) and thus: 

lim C< = {0} . (41) 

/3-s-O 


Another interesting observable is the magnetisation of 
the system for a given energy level Mm = — The 

magnetisation Mm can be easily derived as: 


Mm 


dp 8 Em 
dB 0 dp 


(42) 


From Eq. 


one has: M = 


( dbJg \ 

\dB a ) 


eh . 
2M 2 c 3 ’ 


and 


thus from Eq. (38) the magnetisation My can finally be 




































Figure 2. Magnetisation Mat = — of the system in units of the Bohr magneton. The color schemes is the same as in Fig. jlj 
and follows the energy levels ordered by increasing numerical values of the energy eigenvalue. 


cast as: 


% = ~hb < 


-2N-1—N 2 


1 — AATp ^1 + -^- Af'j 

2N - ±—N 2 

p* 


'l + ANp 1 - — TV 

p* 


p < 0 


p > 0 


(43) 


having introduced the Bohr magneton ps = 2 §jy- 

In Fig. [2]we plot the magnetisation My in units of the 
Bohr magneton as a function of p. It is evident that levels 
from the class C< are characterised by a negative mag¬ 
netisation when p < 0 and by a positive magnetisation 
when p > 0. Exactly the contrary applies to the class C> 
so that in the limit /? —> 0 one recovers the known result 
from ordinary quantum mechanics that the magnetisa¬ 
tion is positive in the left chiral phase (p < 0) while it is 
positive definite in the right chiral phase (p > 0). The 
differences introduced by the new class of states in the 
Anti-Snyder model may have an impact in the termo- 
dynamical properties of such systems and therefore may 
offer a way to access or at least constrain the maximum 
momentum parameter /3. 

Finally we briefly discuss how the results presented in 
this work may have some relevance in the study of new 
materials like graphene [2TIH25] . silicene pTH - lM] and ger- 


manene |291130j in connection with possible experimental 
investigations of the AntiSnyder model. Charge carriers 
in these materials are known to be described by an effec¬ 
tive (2+l)-dimensional Dirac Equation. In reference {37j 
it was shown that both in graphene and silicene the quan¬ 
tum phase transitions induced by the non commutativity 
will affect differently the two inequivalent Dirac points I\ 
and K'. Indeed the quantum phase transitions at K and 
K' arise for different critical magnetic fields. A similar 
behaviour is worth investigating here in the case of the 
AntiSnyder model, although one would have to cope with 
a multitude of critical fields as discussed above. 

In view of the fact that in silicene and germanene the 
charge carriers are massive the exact solutions derived 
in this work can be applied directly while in the case of 
graphene the charge carriers are massless and the same 
solutions would need to be computed explicitly in the 
massless limit, M —► 0. 


VII. CONCLUSIONS 


In this paper we have obtained exact solutions of the 
Dirac oscillator in the presence of a homogeneous mag¬ 
netic held in the Anti-Snyder space. As mentioned before 
the same problem was examined in Snyder space (GUP 
with a minimal length) and it was shown that there were 
a number of quantum phase transitions [BUj accumulat¬ 
ing towards the p = 0 critical point. Unlike the prob- 
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lem in ref. | |35j in the present case the spectrum consists 
of a finite part as well as a continuum (scattering) part 
which will be discussed elsewhere. Similarly to what has 
been found for the same system in the GUP model with 
a minimal length we find that the energy levels can be 
classified in two classes. The levels in one of these classes 
are such that they disappear in the vanishing /3 limit 
(vanishing maximal momentum). Each level in this class 
comes in alternatively in the spectrum (k = 2,4, 6 - - - 
when p < 0 and k = 1,3, 5 • ■ ■ when p > 0). They are all 
characterised by discontinuities which are interpreted as 
quantum phase transitions. The levels in the other class 
instead when (3 —>• 0 reduce to the known spectrum of 
the system in ordinary quantum mechanics, with a single 
critical point at p = 0. 

The magnetisation of the system is an observable that 
in principle has the potential to disentangle the states 
belonging to the class C< from those of the class C>. In¬ 
deed we have found that the two classes are characterised 
by a magnetisation M which is either positive definite or 
negative definite depending on p being positive or nega¬ 


tive. 

It is our belief that the results derived in this work 
may prove useful for further studies of the Anti-Snyder 
model not only from the theoretical side but also from 
the experimental point of view given the recent progress 
in the physics of two dimensional systems involving mate¬ 
rials like silicene, germanene and graphene whose charge 
carriers are effectively described by a (2+l)-dimensional 
Dirac equation. 
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